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INTRODUCTION 


Increasingly, liquid hydrazine (NoH^) is coming to be considered as a 
convenient source of hydrogen rather than just as a rocket fuel. For such 
purposes, the hydrogen is usually obtained by passing the hydrazine through 
a heated catalytic bed. One convenient measure of the effectiveness of a 
catalytic decomposition device as a whole is to compare the quantity of 
hydrogen produced with the equilibrium concentration of the gaseous species 
NH^, N^, and HL which would exist at the temperature and pressure 
found in various parts of the device. Since such data do not seem to be 
readily available in the literature, calculations of the concentrations 
were carried out and are reported here. 

The report begins with a review of the pertinant literature, citing 
experimental results or curve fits to the data where available. The basis 
for the computer program is then described. Following presentation of the 
results in both tabular and graphical forms is a comparison between the 
computed equilibrium concentrations and available experimental data. 






SURVEY OF THE LITERATURE 


This section provides a brief survey of the available literature on the 
thermodynamic properties of hydrazine and its decomposition products as 
well as the decomposition reaction itself. 

Several surveys of the thermodynamic properties of hydrogen, nitrogen, 
ammonia, and hydrazine are available in the literature. The first of these, 
by J. B. Evans and V. J. DeCarlo (1965) was intended as a bibliography for 
investigations into the chemistry of "comet ary phenomena". A second biblio- 
graphy of thermo-chemical data was published by M. K. Buresh, M. L. Reilly, 
G. T. Furukawa, and G. T. Armstrong (1965) of the National Bureau of Stand- 
ards. This particular report listed the low-range temperature range data 
for the heat-capacities, enthalpy, and entropy of ammonia, carbon dioxide, 
and water. Both, reports consisted of references to other investigations. 
Properties of gaseous nitrogen, hydrogen, ammonia, and hydrazine have been 
tabulated and are presented graphically in the N.A.S.A. Design Guide for 
Pressurized Gas Systems, Volume II (1966). 

The data for hydrogen and nitrogen are well known and will therefore 
not be discussed further here. 

J. L. Haws and D. G. Harden (1965) presented empirical equations for 
the thermodynamic properties of hydrazine, including vapor pressure, sat- 
urated liquid density, and the heat capacities. From the 'triple point' 
to the critical point, the vapor pressure is represented by 


(I) InP = 24.24 - (18,1 84.9/T) + 0.47629(Jln T) - 0.003863T + 


(II 15.43(1190.08 - T ) / 1 190. 088 ) &n ( I 1 90 . 8 - T) 

for pressure in psia units and temperature in °R units. The error is 
relatively small. For the heat capacities, these authors give the follow- 
ing relations: 


(2) C p = 0.137858 + 0.052715 x 10 3 T - 0.(19907 x I 0~ 6 T 2 

( 3 ) C v = 0.075903 + 0.52715 x I 0~ 3 T - 0.119907 x lo" 3 ! 3 
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An equation of state for gaseous hydrazine was presented in the form 

-xT/T c „ - xT/T c 

c A, + B^T + C,e c 


(4) P = RT/CV - b) + 


^2 ^ 2 e 


(V - b) 


(V - b)' 


B 5 T 


(V - b) 


(V - b) 


where R 


1 0.73 1 5/mole wt 


b = 0.0299451 ft 3 / I brn 


x = 8.0000 

T = I 1 75.69°R 
c 

A 2 = -32.21898 

B - 0.00294549 

2 

C 2 = -1171.2876 

A 3 = 0.869588 

B 3 = 0. 16838487 x I O -3 

C 3 = 46.250721 

A 4 = -0.202100 x 10” 1 

B 5 = 0.12899897 x JO'" 6 

The equations were used in computing property values for a temperature 
range from 32 °F to I400°F, and for pressures from 1.0 psia to 8000 psia. It 
should be noted that the equations were curve-fitted to match experimental 
data provided by other investigators. 


M. S. Jhon, J. Grosh, and H. Eyring (1967), using the significant 
structure theory of liquids, computed the thermodynamic properties of liquid 
hydrazine (and other chemicals). The computational results were claimed to 
be in good agreement with experimentally-determined property data. These 
properties were I) the vapor pressure, 2) the molar volume, 3) the entropy, 

4) the critical temperature and pressure, 5) the constant volume and constant 
pressure heat capacities, 6) the coefficient of thermal expansion, and 7) the 
compressibility. The surface tension and the dielectric constant were com- 
puted by other methods which are related to the significant structure theory 
of liquids. 
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The "significant structure" theory for liquids is actually a model 
approach for describing the mechanical and thermodynamic properties of 
liquids. 

Let V be the mole volume of the liquid state, V s be the mole volume of 
the solid state at the melting point, and V Q be the mole volume of the 
"sol id I ike structure" in its molten state. The "excess vol ume" of the 
liquid is defined as V - V s . This volume is due to expansion of the chem- 
ical substance upon melting and entering the liquid phase. The liquid is 
considered to have molecular "vacancies" in its interior, and the total 
volume of these vacancies is the "excess volume". This volume permits 
molecules in the liquid state to have translational degrees of freedom. 

If a molecule is completely surrounded by other molecules, then it 
can only oscillate. The oscillation frequency will be nearly the same as 
if the molecule were in a solid state. If the molecule is adjacent to a 
"vacancy", then the molecule will have translational degrees of freedom, 
and will behave as if it were in a gaseous phase. 

The partition function used in the theory is formed from the multi- 
plication of the partition functions for each type of structure. A 
factorial term is included to account for the lack of distinction between 
the 'gas* molecules. 

The liquid state partition function was written in the form 


NV, 


N(V - V Q ) 


( 5 ) 


f = 


[^solid ^degj £fgasj 


NCV - V 0 ) | 


where f^ ^ are the partition functions for solid and gas states, and the 
molecular 'degeneracy of position. 
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The Helmholtz free energy function was expressed as 


(6) 


A = -tfT J2n (f). 


Using the above equation, the vapor pressure was computed in the following 
manner. For a constant temperature, the Helmholtz free energy was plotted 
as a function of volume. Through points common to the liquid and vapor 
states, a tangent was drawn, whose slope gives the vapor pressure. The 
values of volume for the two tangency points give the values for the liquid 
and vapor molar volumes. 

The entropy of vaporization was computed from the relation 




”3A” 


r 3 -i 

(7) 

S = - 

— 

= 

— (kT £n f ) 



_3T_ 

V 

3T 


and the vaporization entropy was given by S - S|, where ’ g T and ’I' indicate 
vapor and liquid states respectively. The constant-volume heat capacity is 
given by 


( 8 ) 


3 

3T 


9 3 &n f 

kT 

3T 


V 


and the constant-pressure heat-capacity is given by 


(9) 


or 


C p = C y t TV 


where 
(10) a 


1 

”3V“ 

1 

< 

Nl 

ro | 

J 

-3^A-j 

V 

_3T_ 

P V 

_3V3T_ 

TV/ 

1 

CM 

> 

CD 

J 


which is the thermal expansion coefficient, and 


(ID 


3V‘ 

3P 



T 


which is the compressibility. The critical temperature, pressure, and volume 
were computed by an iteration method using the conditions 
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( 12 ) 


T 


0 


(13) 


(14) 


3P 

3V 


2 

3 P 




= 0 


rsAn 


p = - 


3V 


The surface tension was computed by an iteration technique which 
accounted for polar liquid orientation effects. The dielectric constant 
was expressed as a function of the index of refraction and the dipole 
moment. 


A comparison of the results of Haws and Harden, with those of Jhon, 
Grosh, and Eyring, for vapor pressure and constant pressure heat capacity 
show discrepancies. Only for the saturated liquid density is there rea- 
sonable agreement. The following table compares the results of the two 
computations for the vapor pressure of I^H^. 


Table I. Comparison of computed vapor pressures. 


T, °K 

P, atm(Haws and Harden) 

P, atmCJhon et_ aj_. ) 

273 

0.003539 

0 

293 

0.0 1 35658 

0,0 1 

313 

0.04322 

j 

0.04 

333 

0.11850 

0.20 

343 

0. 18696 

0.27 


No comparison is possible for temperatures above 400 °K. 

A comparison of the constant-pressure heat capacities gives the 
f o 1 1 ow i ng : 
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Table 2. Comparison of computed constant-pressure heat capacities. 


T, °K 

c p , cal/g-mole-°K 
(Haws and Harden) 

c D , cal/g-mole-°K 
(Jhon et al. ) 

273 

1 J .7792 

20.84 

298 

12.3854 

21.27 

300 

12.4064 

21.35 

400 

14.5728 

j 25.575 


The comparison of results for both vapor pressure and constant-pressure 
heat-capacities show that Jhon, Grosh, and Eyring consistently obtained high- 
er values at elevated temperatures than did Haws and Harden. 

By contrast, the results obtained in each study for the saturated liquid 
density were in good agreement. 


Table 3. Comparison of computed specific volumes. 


T, °K 

V, cc/mole (H&H) 

V, cc/mole (J ,G,&E) 

273 

31 .261 

31 .327 

293 

31 .7269 

31 .67 

308 

32.1294 

31 .984 

400 

35.86127 

34.016 % 387.3 (bp) 


If the results of Haws and Harden, and those of Jhon, Grosh, and Eyring 
are compared with the properties tabulated in the N.A.S.A. Design Guide, the 
results of Haws and Harden are found to compare favorably with those listed 
in the Design Guide. For the purposes of the present study, then, the 
equations of Haws and Harden will be used to provide standards against 
which to compare the results of numerical computations made using the N.A.S.A. 
equilibrium composition analysis computer program. These computations are 
described in detail in subsequent sections of the present work. For the 
moment it is sufficient to note that in these computations each specie is 
assumed to be a perfect gas. 
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Of the major products of hydrazine decomposition, the only one which 
can liquify in the range of temperatures and pressures of present interest 
is ammonia. Since a gas is usually quite imperfect thermally near saturation 
its exact thermodynamic properties under these conditions are of considerable 
interest to one wishing to refine equilibrium composition calculations. 

The thermodynamic properties of ammonia have been tabulated and graph- 
ically presented in the N.A.S.A. Design Guide (1966), Specific properties 
of interest, which are listed, are I) the vapor pressure of liquid ammonia, 
and 2) the density of saturated ammonia vapor, 3) the density of liquid ammo- 
nia as a function of both pressure and temperature, 4) the heat capacities, 
and 5) the heat of vaporization. 

The heat of formation of gaseous ammonia is -11.02 kcal/mole (Penner, 
1968). Equations of state which may be used for ammonia are I) the Berthe lot 
equation, 2) van der Waals, 3) Beattie-Bri dgeman, and 4) the Dieterici equa- 
tion. The equations given below are taken from Penner (1968). The Berthelot 
equation is 

RT 9RT 

(15) V = — + — 

P c 1 28 Pc 

where p c and Tc are respectively the critical pressure and temperature. The 
gas constant R for ammonia is obtainable, for any set of units, from the 
N.A.S.A. Design Guide. 

The van der Waals equation would give 



(16) P + v * b j = RT 

where a = 3p c V c 2 and b = V c /3, 

The Beattie=Bri dgeman equation for gaseous ammonia is 



where A 0 = 2.3930 

B q = 0.03415 

a = 0.17031 

b = 0.019112 

c = 476.87 x I0~ 4 
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for units in atmospheres, liters per mo|e, and °K. 


The Dieterici equation, which is very accurate near critical points, is 


(18) 


Pe 


VRT 


(V - b) 


where a - 7.39p c V c and b = V /2. 


RT 


The values for the constants in these equations are from Hi rschfelder, 
Curtiss, and Bird (1954). 


The remainder of this review is concerned with papers treating the 
decomposition of hydrazine. 


K. W. Michel and H. Gg. Wagner (1965) studied the thermal decomposi- 
tion of diluted with argon and/or helium, behind shock waves, in the 

temperature range MOO °K to 1600 °K, for pressures from 26.65 atm to 80.0 
atm. For temperatures between MOO °K and 1200 °K, decomposition did not 
start until an 'induction period' had passed. For temperatures greater than 
1300 °K, the decomposition reaction begins by following a first-order reac- 
tion-rate law, and then changes to successive reactions of varying order. 

A reaction of order one follows a rate law such as 

dc 

(19) - — = kic 

dt 


where c = specie concentration, k| = the rate constant, and t = time. 
Reactions may have an order greater than one. An example of such a reaction 
is 


(20) 


A + 


X + 


Y . 


Let c a be the concentration of specie A and eg be the concentration of specie 
B. Let the superscript o indicate the value of a quantity at time t = 0. 

The reaction rate equations are 


( 21 ) 


Def i ne x 
(22) 



dc A 

dCj 


— — 

dt 

dt 


k 2 c A c B ‘ 



C B 


dx 

— = k 2 (c A ° - x) (c B ° - x) - 
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Cg° , then 


If, as a special 

case. 

r* c 

C A 


dx 


(23) 

— 

== 


dt 


(24) 


d° A 



dt 


o - 


2 A 


- k 2 c A . 


Many reactions are not as simple as those given above, and the reaction 
order may be non-integer in value. 

The results of Michel and Wagner's investigation confirmed the "straight- 
chain" model for the decomposition reaction which had been proposed by G. K. 
Adams and G. W. Stock (1953). Reaction models for the hydrazine break-down 
are discussed later. At higher temperatures, the reaction rates were found 
to be partial I y^dependent on the total gas density. For a density of 7.5 
x I0 “ 5 moles/cm , (26.65 atm), 

(25) log k = 12.8 - ( (52,000 cal/mole)/2.3RT ) 

—5 3 

and for a density of 2.5 x 10 moles/cm , (80.0 atm) 


(26) log k = 12.0 - ( (48,000 cal/mole)/2.3RT ), 

The main objective of the investigation was to determine the effects of 
small amounts of oxygen on the decomposition rate. It was found that small 
amounts, (i.e. 0.2$ of the mixture) do not affect the reaction. For temper- 
atures above 1400 °K, heavy concentrations of oxygen, (e.g. =. 10) 

still do not affect the decomposition reaction. 

The amount of ammonia formed during the decomposition was also deter- 
mi n 
NLH 
to 

R. Roback (1965) made an investigation into the thermodynamic properties 
and equilibrium compositions of the products of dissociation of hydrogen, 
ammonia, and other materials, such as graphite, molybdenum, tungsten, and the 
carbides, nitrides, and oxides of titanium and zirconium. The temperature 
range was from 1800 °R (1000 °K) to 10,000 °R (5550 °K) for a pressure of 
1000 atm. The materials considered were for use in controlling the radiation 
heat transfer to the moderator walls of gaseous nucl ear-engi ned missies. The 
computations were carried out by assuming that all species obeyed the ideal 
gas law. Viscosities and thermal conductivities were obtained using data 
valid at one atm pressure, and a pressure correction determined from the 
Enskog theory for dense gases. Atomic and molecular species were considered 


ed . It was found that at 1100 K the ratio of NH 3 formed per amount of 


decomposed is equal to one, at 1600 °K equal to 0.5, and at 2000 
Eero. The radical NH appeared for temperatures above 1400 °K. 


K equal 



in the computation, but ionic species were considered to be of negligible 
importance below a temperature of 10,000 °R (5550 °K). The calculations 
indicated that hydrogen dissociation is important for temperatures above 
7000 °R (3880 °K). By contrast, concentrations of ammonia are almost 
completely decomposed into H 2 and N 2 at 2000 °R (1110 °K) . Densities and 
enthalpies were computed for hydrogen, ammonia, and other gases, and pre- 
sented as functions of temperature. 

I. J. Eberstein and I. Glassman (1965) experimentally studied the 
decomposition of low concentrations of hydrazine for a temperature range 
of 750 °K to 1000 °K. The apparatus used was an adiabatic flow reactor 
with a driving pressure of 136 atm. Both the rate constants and reaction 
stoichiometry were determined. For hydrazine break-down inside a 3-inch 
duct, the reaction rate constant was determined to be 

(27) k = I0 i0 * 33 exp(-36,I70/RT) sec -1 . 

The over-all reaction, of an order close to one, was given as 


(28) N 2 H 4 ► 0.9NH 3 t 0.5N 2 + 0.6H 2 . 

The stoichiometry of the reaction was determined by experiment. 

Hydrazine-water mixture decomposition was also studied. For mixtures 
where the water was more then one percent of the total volume, the rate of 
decomposition was decreased by a factor of ten. The reaction rate for these 
mixtures was found to be independent of the amount of water present in the 
mixture, provided that the amount of water was more than one percent of the 
total volume. It was concluded that water must suppress one of the decom- 
position reaction steps. 

Along with hydrazine, both unsymmetrical dimethyl hydrazine and mono- 
methyl hydrazine were considered. Of the three, anhydrous hydrazine, N 2 H 4 , 
had the slowest decomposition rate. 

E. J. Bair et aji_. (1966) used the explosive decomposition of hydrazine 
as part of their experimental program in the development of a high-resolution 
spectrometer and the spectroscopic procedures in examining explosive reac- 
tions. The reaction was described as a series of "free radical chain pro- 
cesses". One of the first products to be formed was NH 2 , but this product 
disappeared as the reaction continued, forming NH. The NH radical in turn 
also disappeared. Both radicals were observed spectrometrical ly. When the 
NH radical vanished, a spectrum of lines in the 0.600p region was found 
which was thought to arise from the formation of N 2 H 3 , termed the "hydrazil 
rad i ca I " . 
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A. S. Sokol ik, V. P. Korpov, and E. S. Semenov (1967) used hydrazine 
decomposition flames in their study of turbulent burning rates and flame 
velocities. The decomposition process was considered to be a "straight-chain 
mechanism." The breaking of the reaction chain, however, is dependent on 
pressure. The reaction order is also dependent on reaction pressure. At 
low pressures, the order is approximated I y 1.0, and at high pressures the 
order is 1.5. It was noted that hydrazine "decomposes extremely readily," 
requiring special care in both the handling and experimental procedures. 

R. F. Sawyer and I. Glassman (1967) have examined the reactions between 
hydrazine and nitrogen dioxide, nitric oxide, oxygen, and oxygen-nitric 
oxide, by using an adiabatic flow reactor. For each reaction, the investi- 
gators determined reaction rates, activation energies for the over-all 
reaction, and the reaction orders. Reacting with oxygen, hydrazine forms 
water and nitrogen according to the reaction 


(29) N 2 H 4 + 0 2 2H 2 0 + N 2 

For a total pressure of 1.0 atm and temperatures between 950 °K and 
1010 °K, the hydrazine/oxygen reaction rate law is expressed as 


(30) d[N 2 H 4 3 = - kCN 2 H 4 ] 

with dt 

(31) k = 1 0 9 * 9 1 exp ( -37 , 200/RT ) sec" 1 

and a reaction order of one. The activation energy was 37.2 kcal/mole. 

The reaction rate is dependent on the concentration (to power one) 
and independent of the concentration of oxygen. Oxidation was shown to be 
possible whether or not the hydrazine was decomposed prior to the reaction. 

The production of water by the reaction was believed to be responsible for 
slowing the rate of decomposition. 

D. I. Maclean and H. Gg. Wagner (1967) studied the decomposition-flame 
structure for both ammonia -oxygen and hydrazine vapor, determining experi- 
mental ly temperature and specie-concentration profiles. The initial tempera- 
tures of the hydrazine vapor were between 70 °C and 90 °C, for pressures of 
0.01842 atm and 0.0342 atm. It should be noted that the only initial sub- 
stance in combustion was the pure hydrazine, with no other chemical species 
present. The reaction, after ignition, proceeds without aid of any oxidizers, 
and is exothermic. For the hydrazine-vapor flame, a maximum temperature of 
1350 °K was found by means of a thermocouple. The investigators believed the 
temperature reading to be too low, as the adiabatic flame temperature should 
be near 1700 °K. 
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It was found by experiment that increasing the flame speed decreased 
the rate of decomposition of the first few millimeters above the burner. 

A major factor in the rate of decomposition is pressure. By doubling the 
pressure, nearly twenty mole percent of the hydrazine was broken down within 
the first millimeter above the burner. 

The concentrations of H 2 , N 2 , and NH 3 , which are the products of the 
reaction, and the rate of the reaction, are not affected by changes in the 
initial temperature of the hydrazine vapor. The maximum temperature is not 
affected by changes in pressure. 

E. F. Logan and J. M. Marche 1 To (1968) decomposed hydrazine by means of 
an alternating-current glow-discharge flow reactor. Using several flow rates 
and electric currents, light emmission spectra were measured from 0.300 p to 
0.600 p. The most intense spectrum was produced by molecular nitrogen. How- 
ever, atomic nitrogen spectra were not observed. Both atomic and molecular 
hydrogen were observed. The emmision bands for ammonia were not found. The 
radicals NH, NH 2 , were found but OH was not. The gaseous hydrazine vapor 
pressure was 1.91 x I0 -2 atm. 

It was found that 65$ of the hydrazine is decomposed very rapidly, with 
the rest of the reactant decomposing at a slower rate. 


The decomposition of hydrazine has also been used in the study of the 
dissociation energy of the N-H and N-N bonds. 

I. P. Fisher (1965a, b), with G. A, Heath, experimentally determined 
the value for the dissociation energy of the N-H bond. The value was found 
by electron- impact measurements, and by using the relation 


(32) 


D(H - N 2 H 3 ) = A(N 2 H 3 + ) - l(N 2 H 3 ) 


where A is the appearance potential of N 2 H 3 + , I is the ionization potential 
of N 2 H 3 , and D is the dissociation energy of the N-H bond. The values which 
were found were 


(33) 

(34) 

(35) 


A(N 2 H 3 + ) 
I (N 2 H 3 ) 
D(H-N 2 H 3 ) 


11.3 - 0.1 ev 
7.88 ev 

78.0 kcal/mole , 


Values for bond strengths, and heats of formation and ionization potentials 
of several chemical species resulting from hydrazine decomposition were also 
determined (Fisher, 1965b). 
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Another study of bond dissociation energies, for the N-N bond, was made 
by H. A. Olschewski, J. Troe, and H. Gg. Wagner (1966). It was shown that 
the decomposition activation energy must be measured for decomposition at 
high pressures. A shock tube of 70 mm diameter was constructed and used in 
their experimental program, with driving pressures up to 200 atm being 
app lied. 

For decomposition at high pressures, and densities greater than 5 x 10“^ 
moles/crn^, the experimental results provided the temperature dependence of 
the reaction rates. The reaction rates were then used in determining the 
uni molecular decomposition activation energy, which was found to be 55 - 2 
kcal/mole. The reaction rate for decomposition at high pressures was deter- 
mined to be 


(36) k = l0 l3 - 6 exp(-55,000/RT) sec ’. 

This rate constant is for the over-al I reaction, and not to be confused with 
the dissociation energy of the N-H bond. 

Two papers were found which discussed the decomposition of hydrazine by 
catalytic means. 

A. S. Keston (1967) carried out a theoretical study of catalytic react- 
ors designed for hydrazine decomposition.* Transient and steady-state models 
were developed for temperature and specie-concentration profiles within the 
reactor, as functions of time and position. Both thermal and catalytic- 
induced decomposition were considered, with the inclusion of the effects of 
heat and mass transfer between the free-flowing gas, and gas contained within 
the catalyst material. The results of the computations were compared with 
data from hydrazine-fueled engine experiments where the Shell 405 catalyst 
was used, with very good agreement between the theoretical and experimental 
results. These models provide an indirect means by which the decomposition 
of hydrazine and ammonia may be computed. 

A summary of the results of the steady-state reaction calculations of 
Keston is given on the following page. The values are for an axial distance 
of 0.25 feet from the inlet end of the catalytic reactor. 


Temperatures varied from 1000 °K to 1330 °K and pressures varied from 7.6 
atm to 71 atm. 
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Table 4. Summary of steady-state reaction computations of Keston 


p 

T 

Mole Flow Rate 

Moh 

3 Fractions 

Fractional 

atm. 

°K 

1 bm/ft^-sec 

h 2 

N 2 

nh 3 

NH3dissoc. 

7.59 

1000 

1.51 

0.55 

(0.59) 

0.31 

0.14 
(0. 10) 

0.62 

14.8 

1020 
( 1000) 

1.52 

0.57 

(0.59) 

0.32 

(0.33) 

0. 10 
(0.07) 

0.71 

32.6 

1054 

(1076) 

3. 12 

0.55 

(0.56) 

0.31 

0.14 

0.64 

66.2 

1 120 

6.29 

0.53 
(0.51 ) 

0.31 

(0.32) 

0.16 

(0.17) 

0.57 

71 .0 

1054 

(1076) 

2.43 

0.57 

(0.54) 

0.32 

0.12 
(0. 15) 

0.69 


The numbers in ( - ) represent experimental data differing from the theo- 
retical computations. 

The mole fractions of H 2 and N? produced by the reaction are nearly 
independent of variations in pressure and temperature. The mole fraction 
of hydrogen produced varied between 0.53 and 0.57, and the mole fraction 
of nitrogen produced varied from 0.31 to 0.32. Ammonia mole-fractions 
varied between 0.10 and 0.16. It was found that in the reaction, mole- 
fractions of ammonia reached a maximum with respect to axial distance in 
the reactor (or equivalently, time), and then declined In value. Between 
57$ and 71$ of the NH 3 produced in the reaction decomposes to nitrogen 
and hydrogen. 

It may be noted that for the maximum mass flow rate, hydrogen production 
was at a maximum, as was the fractional NH 3 dissociation. However, the tem- 
perature of the reaction reached a maximum for lower mass-flow rates. 

At either 14.8 atm or 71.0 atm, there was maximum output of hydrogen, 
chiefly due to the maximum dissociation of ammonia. 

A second study of catalyst- induced hydrazine break-down was made by 
C. F. Sayer (1970) who also used the Shell 405 catalyst. The purpose of 
the study was to examine the chemical kinetics and mechanisms involved 
in the break-down of liquid hydrazine in the presence of the catalyst. 
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Using the Shell 405 catalyst, the break-down resy|ted in a gas mixture 
which was reported to be 99.4$ N 2 and 0.6$ N 2 H 4 . The vanishing of the 
hydrogen was attributed to " the conversion of the hydrazine molecule into 
the nydrazinium ion." If a rhodium catalyst was used instead of the Shell 
405 catalyst, then the resulting gas mixture was nearly 50$ N 2 and 50$ F^. 
It may be noted that a main component of the Shell 405 catalyst is iridium. 
The reaction with the Shell 405 catalyst was expressed as 


(37) 


3N 2 H 4 



+ 4NH, . 

-J 


For the rhodium catalyst, the reaction was written as 


(38) 


2N 2 H 4 


2NH 3 . 


For the production of hydrogen from hydrazine, the rhodium catalyst is 
obviously to be preferred. 

Two references attest to the volatility of hydrazine. G. P. Sutton 
(1963) noted that hydrazine may spontaneously ignite with either nitric 
acid or hydrogen peroxide. It was also noted that hydrazine vapor may 
form explosive mixtures with air. 


Sutton very briefly discussed the reaction of hydrazine with catal- 
ysts. The products of the decomposition reaction, the exact composition 
of the product mixture, and the pressure will depend on the type of catal- 
yst used. 


If impurities are present in liquid hydrazine, or in contact with 
hydrazine vapor, then the hydrazine will readily decompose exothermic- 
ally. 


The reaction was stopped when the total pressure in the experimental 
apparatus reached 1.025 atm. 
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E. J. Bowen and A. W. Birley (1951) found by experiment that for very 
low pressures, and temperatures of the order of 100 °C, that mixtures of 
hydrogen peroxide and hydrazine vapors may co-exist for a brief period of 
time without reaction. However, decomposition will then take place. 

A comparison among the results for hydrazine decomposition shows fair 
agreement between the investigations of Michel and Wagner (1965), Eberstein 
and Glassman (1965), Sawyer and Glassman (1966), and Maclean and Wagner 
(1966). Without the presence of a catalyst, the stoichiometry of the 
decomposition reaction is generally 


(39) 


N 2 H 4 



iH 2 



* 


By inducing decomposition with a glow-discharge at a pressure of 1.91 
x I0” 2 atm, Logan and Marchel lo (1968) obtained a gas mixture consisting of 
46% H 2 , 28% N 2 , and 19% NH^ for 85% hydrazine break-down. 

C. F. Sayer (1970) found that catalyst-induced reactions give results 
different from the previous investigations mentioned above. Using the Shell 
405 catalyst, one mole of hydrazine results in the formation of 1.33 moles 
of ammonia. The final gas composition is however 99.4% N 2 , and 0.6% N 2 H 4 . 
If a rhodium catalyst is used, one mole of hydrazine decomposes into one 
mole of ammonia, with a final gas composition of 50% N 2 , and 50% H 2 . The 
chemical reactions were stopped when total system pressures exceeded 1.0 
atm. 


Reaction pressures varied considerably. The pressures in the experiments 
were 


Michel and Wagner 
Eberstein and Glassman 
Sawyer and Glassman 
Maclean and Wagner 


26.26 atm to 80.0 atm 
136 atm 
1 atm 

0.0184 atm to 0.0342 atm 
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A, S. Keston (1967), by theoretical means, found that for elevated 
pressures, the results of decomposition of hydrazine are mostly Ho and N 2 . 

As an example, for a pressure of 7.6 atm, and a temperature of 1000 °K the 
gas mixture consisted of approximately 54$ H 2 . 31$ N 2 , and 15$ NH,. The 
maximum temperature was 1120 °K, and the maximum pressure was 71.0 atm. 

For relatively low pressures, and temperatures of the order of 1000 °K, 
the present numerical investigation shows that the products of the decom- 
position reaction are almost entirely hydrogen and nitrogen, with twice as 
many moles of hydrogen being produced as nitrogen. 

In comparing these results with those of previous investigations, the 
agreement is not unreasonable. In previous studies, for very low pressures, 
one mole of hydrazine decomposed to approximately one mole of ammonia, and a 
half-mole each of hydrogen and nitrogen. If a rhodium catalyst is used, the 
products of the reaction are equal mole numbers of hydrogen and nitrogen, with 
no ammonia left at the end of the reaction. 

For elevated pressures (Keston, 1967), hydrogen is the major product 
(54$) with lesser amounts of nitrogen (31$) and ammonia (15$). It may be 
seen that as pressure increases, for temperatures in the neighborhood of 
1000 °K, the ammonia production declines, and hydrogen and nitrogen are 
produced in amounts which appear to approach a ratio of two to one. The 
present computational prediction is that for pressures in the neighborhood 
of 7.6 atmospheres and temperatures of about 1000 °K, the mole-fractions 
of the decomposition products are 34$ nitrogen and 66$ hydrogen. 


Several reaction mechanisms have been postulated for the break-down of 
hydraz i ne. 


Michel and Wagner (1964) have noted that the first two reaction steps 
have been theoretically determined and experimentally verified. These steps 
are 

(40) N 2 H 4 — ► 2NH 2 

(41) NH 2 + N 2 H 4 — ► N 2 H 3 + NH 3 . 

For low temperatures, the next steps would then be 

(42) N 2 H 3 ► H + HNNH 

fol lowed by 

(43) H + N 2 H 4 NH 3 + NH 2 
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or 


(44) 

H + N H ► 

H + N 9 H, 


2 4 

2 2 3 

(45) 

NH 2 t HNNH — - 

NH 3 + N 2 + H 

(46) 

2NH 2 — — ► NH 3 

t NH 

(47) 

M + NH 2 + H 

► NH 3 t M , 


At temperatures of 1600 °K, the step 

(48) 2NH — HNNH + H 2 

preceeds reaction (6). For temperatures in the neighborhood of 2000 °K, step 
(6) could be replaced by the reaction 


(49) HNNH NNH + H - 

Eberstein and Glassman (1965) considered two possibile reaction mechan- 
isms. The simple mechanism was similar to the postulated mechanism of Michel 
and Wagner, except for the last step. Michel and Wagner assumed the reaction 


(50) NH 2 + HNNH — — ► NH + N 2 + H 

whereas Eberstein and Glassman consider the reaction to be 


(51) N 2 H 3 + NH 2 — ► NH 3 + N 2 + H 2 . 

The assumed mechanism includes branching reactions. By numerical 
techniques, the investigators found that the reaction has four parts. In 
the first part, the formation of free radicals takes place. The reaction 
steps which follow can be characterized by "steady-state" rate constants. 
The third part of the reaction was characterized by decreasing first-order 
rate constants. In the last part, the concentrations of free-rad i ca I s and 
the over-all reaction rate rapidly declined. When compared with experi- 
mentally-measured stoichiometry and reaction-rates, the reaction mechanism 
was found to be "in good quantitive agreement." 

The reaction mechanism which was postulated is tabulated as follows. 
Initiation: 


(52) 


n 2 h 4 


t X 


2NH, 


+ X . 
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Propagation: 


(53) 

(54) 

(55) 

Branch i ng: 

(56) 

(57) 

Termi nation: 


(58) 

nh 2 + 

n 2 h 3 - 

— ► NH 3 + n 2 + 

(59) 

N„H^ + 

N 0 H, - 

— ► 2NH 3 + N 2 


2 3 

2 3 


(60) 

n 2 h 3 + 

H ■ 

► N 2 + 2H 2 

(61 ) 

NH 2 + 

nh 2 1 

- n 2 h 4 . 


Sawyer and Glassman (1966) have tabulated the assumed mechanisms for 
reactions of hydrazine with nitrogen dioxide, nitric oxide, and oxygen. In 
the reactions given, the authors assumed that branching occurs before propa- 
gation. 

Sokol ik, Korpov, and Semenov (1967) have noted that the "chain" mechan- 
ism of Adams and Stocks is usually accepted. The steps for this mechanism 
are 


N 2 H 3 + X ► NH + NH 2 + X 

N 2 H 4 + NH — -*■ NH 2 + N 2 H 3 , 


n 2 h 4 + nh 2 — ► nh 3 + n 2 h 3 

N 2 H 3 + X ** N 2 + H 2 + H + X 

N 2 H 4 + H NH 3 + NH 2 . 


(62) N 2 H 4 ► 2NH 2 . 

Continuation: (propagation) 

(63) N 2 H 4 + NH 2 ► NH 3 + N 2 H 3 

(64) N 2 H 3 ► N 2 + H 2 t H 

(65) H + N 2 H 4 . ► NH 3 t NH 2 . 
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The chain breaking is dependent on the reaction pressure. For low 
pressures, the branching and termination steps are 


(66) H + N 2 H 3 •*- N 2 + 2H 2 

(67) NH 2 + N 2 H 3 ► NH 3 + N 2 + H 2 . 

For high pressure reactions, the branching and terminal steps are 

(68) H + NH 2 + X ► NH 3 + X 

(69) H + H + X ► H 2 + X . 

Logan and Marche I io (1968), discussing the glow-discharge break-down of 
hydrazine, gave the following reaction steps as the mechanism of electron- 
induced decomposition. The first step is caused by the electron interaction 

(70) N 9 H, + e* <*- N_H * + e 

^ Z 4 

The (*) indicates an activation state. The authors' spectrometer measure- 
ments indicated that N 2 H 4 * breaks down in two simultaneous reactions. 

(71) NJH * ► N H + H 

2 4 2 3 

(72) N 2 H 4 * ► NH 3 + NH 

The reactions which then fol low are 

(73) H + N 2 H 3 . H 2 + N 2 H 2 

(74) NH + N 2 H 4 NH 3 + N 2 H 2 . 

The diazene (or diimide) decomposes as 

(75) N 2 H 2 ► N 2 + H 2 . 

In all of the assumed reaction mechanisms, N H always participates in 
a later reaction step, as well as in the initial step. There is disagreement 
as to the exact steps of the reaction after the initial and propagation reac- 
tions. 

Sayer (1970) formulated a hydrazine break-down mechanism involving the 
use of catalysts, to explain the results of his experiments. It was postu- 
lated that hydrazine molecules adsorb on the catalyst surface, and react with 
other hydrazine molecules which are "in solution" to form two "ci s-hydrogen" 
atoms. Additional products of this reaction are two molecules of ammonia, 
with cis- diimide adsorbed on the catalyst surface. Two adjacent molecules of 
diimide then react to form nitrogen and hydrazine. The nitrogen desorbs from 
the catalyst surface, and is replaced by another hydrazine molecule. 
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H. S. Maddix (1965) discussed the break-down of ammonia in his study of 
gaseous discharge phenomena. The ammonia decomposition was induced by an 
electrical discharge in a tube. The observed reaction was 


(76) 2NH 3 N 2 + 3H 2 . 

As the amount of ammonia decreased, the nitrogen concentrations steadily 
increased. The hydrogen concentrations first rapidly increased, the de- 
creased, and then increased again. The changes in the amount of hydrogen 
were due to chemical adsorption of the hydrogen with the tube walls. 

A second discussion of the decomposition of ammonia was provided by 
Roback (1965). It was noted that above 2000 °R, ammonia is almost entirely 
dissociated, the products being H 2 and No. Above 4000 °R, H, NH, and NH 2 
concentrations become appreciable. For temperatures above 7000 °R, N 2 
begins to decompose to single atoms of nitrogen. 

Michel and Wagner (1965) included a study of the pyrolysis of NH , which 
showed that ammonia decomposition is a first-order reaction. For tempera- 
tures between 2100 °K and 2900 °K, and with a ' carrier’ gas of argon, the 
reaction rate was 

(77) k = IO l5,64 exp[-(79.5 - 2.5 kcal/mole)/RT] cm 3 /mole-sec . 

For temperatures greater than 1400 °K, NH radicals will be present in the 
decompos i t i on-products mixture . 

As has been shown, values of thermodynamic properties of hydrogen, nitro- 
gen, ammonia, and hydrazine are readily available. The properties of hydra- 
zine have been determined from both expert mental studies and theoretical com- 
putations. Empirical equations for the properties of hydrazine have been 
presented, and equations of state for both gaseous hydrazine and ammonia have 
been listed. 

In the forego! ng, previous studies of hydrazine decomposition have been 
summarized and discussed. The studies have shown that in general one mole 
of hydrazine decomposes into one mole of ammonia, one-half mole of hydrogen, 
and one-half mole of nitrogen. The exact stoichiometry of the reactions 
depends on both the temperatures and pressure of the reaction. The use of 
a catalyst can lead to increased production of both hydrogen and nitrogen. 
However, if the incorrect catalyst is used, very little hydrogen will result 
from the hydrazine decomposition. 

The results of the previous investigations and present numerical study 
have, where possible, been compared. 
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As has been shown, the decomposition mechanism is quite complex, and 
there is still disagreement among researchers studying the hydrazine break- 
down as to the exact steps, in the reaction. Several of the investigations 
have included studies of ammonia decomposition, and these have been sum- 
marized. 

It is clear that the properties of hydrazine are not as well determined 
as one would wish. Although previous studies are not in exact agreement, 
there is qualitative agreement among them. Much more experimental work, 
however, will have to be performed to determine the exact kinetics of the 
decomposition reaction mechanism. 
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THE COMPUTATIONAL METHOD 


The computer program used in the present calculation of hydrazine decom- 
position was written by S. Gordon and B. J, McBride (1971), and listed in 
N.A.S.A. SP-273. The program was designed to compute the equilibrium com- 
position of reacting chemical mixtures for given thermodynamic states, i.e. 
specified (T,P), <H,P), (S,P), (T, V ) , (U,V), or (S,V), where: 


T = temperature, °K 
P = pressure, newton/rneter^ 

V = specific volume, mete rV kg 
H = enthalpy, joules 
S = entropy, joules/°K 
U = internal energy, joules. 

The chemical species may be either in a gaseous phase and/or in a condensed 
phase. 

The computer program also has the capability of computing rocket 
performance parameters, incident and reflected shock parameters, and Chapman- 
Jouguet detonation properties. 

The method used for the present set of equilibrium composition calcu- 
lations is based on the concept that reacting chemical systems are in equi- 
librium when the free -energy of the system is at an extremum value. The 
specified thermodynamic variables in this report are always pressure and 
temperature. The appropriate thermodynamic function for these variables 
is the Gibbs free energy function. As shown by J. B. Cal ten (I960), this 
function will possess a minimum. The minimization of the Gibbs free energy 
is constrained by additional physical conditions. For chemical reaction 
systems, one appropriate condition is that the mass involved in the reactions 
be conserved. When a given function possesses an extremum value, which is to 
be determined under one or more restraining conditions, the most suitable 
technique for determining this is the method of Lagrangian multipliers. 

The alternative method of calculating chemical equilibrium composi- 
tion is by using sets of simultaneous equilibrium equations in terms of 
equilibrium constants. The differences between the two methods, with 
their advantages and disadvantages, has been discussed by F. J. Zeleznik 
and S. Gordon (i960, 1966, 1968), and S. Gordon (1970). The invest iva- 
tions and results reported in the above-cited papers formed the basis 
for the computational method and program outlined by S. Gordon and B. J. 
McBride (1971 ). 
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Gordon (I960, 1970) has shown that if ’generalized* methods of solution 
are used to solve equilibrium problems expressed in either free-energy terms, 
or by equilibrium constant equations, the resulting sets of equations will 
reduce to the same number of equations to be solved. However, the ’constant' 
method has disadvantages not present in the free-energy method. The first 
of these, as noted by Gordon (1970), is "more bookkeeping". In using the 
’constant’ approach, the components (amounts) of the reactants and products 
must be selected, and the possible reactions written in terms of these com- 
ponents. Furthermore, equilibrium constants must be computed. It may be 
noted that such constants are not always available in the literature for 
high pressures and/or temperatures. The free-energy method does not require 
any part of the above-mentioned procedures. 

If a chemical species is present in only trace amounts, i.e. "small 
components," then problems will occur in numerically solving the equilibrium 
constant equations. In the free-energy method, numerical difficulties are 
easily circumvented. 

When one or more species is present in the form of a condensed phase 
(liquid and/or solid), the composition problem is very difficult to solve 
using the ’constant’ approach, but is much easier to solve by the free-energy 
technique. 

If non-ideal equations of state are to be used, then the 'constant' 
method will be very difficult to use, as the generalized solution methods 
are not very amenable to such modifications. 

Partition functions may also be used in solving equilibrium problems. 

As shown by Lee, Sears, and Turcotte (1963), the equilibrium constant for 
a given chemical reaction may be expressed in terms of the partition functions 
of the reacting species. As an example, consider the reaction 


( I ) X t Y = XY . 

The equilibrium constant for the reaction is 

k = Jzi. 

Z X Z Y 

Z x = the partition function for specie X 

Zy = the partition function for specie Y 

Z^y = the partition function for specie XY , 


( 2 ) 

where 
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The Gibbs free energy function can also be written in terms of the 
partition functions: 


(3) 


-k'T£n 



OX 

c 

— 

n° 

Y 


fx 

o 

Z Y 

o 

Z XY 

o 

n x 

-J 


n Y_ 


n XY 



where k’= the Boltzmann constant, and n^, n°, and n^y are the number of 
molecules of species X, Y, and XY at equilibrium. 


Using the partition functions in the calculation of equilibrium chemical 
compositions involves computing equilibrium constants, a step which is not 
necessary in the f ree-energy method. 


in the N.A.S.A. computer program, it is assumed that all gaseous 
chemical species will behave as ideal gases. For the high temperatures and 
pressures considered in the present work, the ideal gas law will be a fair 
approximation to the real behavior of the reacting system. The equation of 
state for ideal gases is written 


P 

(4) - = NRT 

P 

m 

where N = £ n ■ 

nj = the number of kg-moles of the j"*"* 1 species. 

m = the number of species in the gas mixture 

n - m = the number of condensed species . 

The following discussion is based in part on the material in the reports 
of Zeleznik and Gordon (1968) and Gordon and McBride (1971). The Gibbs free- 
energy function may be defined as 


(5) G-f - Gf (T,P,nj ) 
or equivalently, 

(6) G f = U - TS + PV 
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and as 


m 

(7) dU - TdS - PdV + 2 U-dn. 

j J J 


then 

( 8 ) 


m 

dGf = -SdT + VdP + v ^j dn j 

j 


where the chemical potential (per kg -mole) is defined as 


(9) 


l- 1 ! 



3n , 


T,P, n 


k / j 


The Gibbs function is homogeneous in n . , 
the property that 


i.e. the Gibbs function has 


(10) Gf (T,P ,An . ) - AG^ (T, P, n . ) « 

As previously mentioned, the Gibbs free energy function is appropriate 
in analyzing systems where the pressure and temperature are of constant value. 
It was also stated that the function will have a minimum value for a system 
in equilibrium. A system is considered to be in equilibrium if it is not 
changing, either mechanical I y , thermally, or chemically. 

To show that the Gibbs function of a system is of minimum value when 
the system is in equilibrium, consider a system in contact with two ’reser- 
voirs’, one at constant temperature, the other at constant pressure. The 
general conditions for equilibrium are that 


(11) dU = 0 

(12) d 2 U = 0 . 

For an equilibrium situation, the total energy of the given system and 
the two ’reservoirs' must be of minimum value. Hence, 


(13) d(U + U p ) = 0 

where r refers to ’reservoir’ conditions. For the constant temperature 
reservo i r. 


(14) 


dU r ^y = -TpdS 
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and for the constant pressure reservoir 


(15) dU p = P r dV. 

Therefore, 

(16) d(U + U ) = dU - T r dS + P dV 

r 1 r 

(17) d(U + U r ) = d(U - T r S + P p V) . 

Since the given system will have a temperature equal to 1 r , and a 
pressure equal to P , 

(18) d(U + U r ) = d(U - TS + PV) 

(19) d(U + U p ) = dG f 

and the condition for equilibrium becomes, for constant temperatures and 
pressures. 


(20) dG f = 0 . 

For chemically reacting systems, the equivalent condition is determined 
by using the Euler relation (Callen, I960) and substituting into equation 
(18). 

m 

(21) U = TS - PV + 2 y j n j 

j = I 

The equilibrium condition then becomes 

m 

(22) d( 2 y.n.) = 0 . 

j = I 

To show that the condition is for a minmum value of the Gibbs function, 
consider the second derivative of the system energy. 

(23) d 2 (U + U p ) = d^ + T p d 2 S + P p d 2 V = d 2 G f 
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and for the reservoirs. 


(24) T_d 2 S = 0 

(25) P r d 2 V * 0 
and from equation ( 12 ), 


(26) d 2 (U + U p ) - d 2 G f > 0 • 

The Gibbs function thus has a minimum value for a system in equilibrium. 

The free-energy per kilogram of gas mixture is, for constant temperature 
and pressure. 


(27) 9 - • 

j = I 

The condition for the mass balance of the reaction system needs to be con- 
sidered. The expression for the mass balance is 


(28) 



£ 


where 



the stoichiometric coefficients (the kg-atoms of 
element i per kg-mole of specie j ) 


the specified number of kg-atoms of element i per 
kg of total reactants . 


Define a function G as 


(29) 


where 


g = g 



V b i - V 


I = 



A 5 = a Lag rang Ian multiplier 
% = the number of chemical elements# 
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Equation (29) represents the Gibbs Function, which is to be minimized, 
with the applied constraint. The condition for equilibrium becomes 


(30) 


n £ 

6G = 0 = 2 Pi + ? Xtar : 

j=l J i-i J 


6nj +2 (b { - bVJdAj . 

i = I 


If the variations in n. and A. are independent, i .e. independent 6n. and 
6Aj , then the resulting equilibrium equations are J 


£ 

(31) Pj + £ ^ i a j j = 0 

i = I 

(32) bj - b° = 0 . 

It should be noted that A contains the condition for phase equilibrium, 
and that the chemical potential of a species is the same for all phases of 
the species at a given temperature and pressure. 

For the ideal gas equation of state, the chemical potential of the 
species may be expressed as 


(33) y , = y? + RT £n 

v J 

j = I , . . . , m 

(34) yj = y? j = m + I , . . . , n 

where yj is the chemical potential in the standard state, and P g ^ m is 
pressure given in terms of atmospheres. r l! 


+ RT £n 


(P a+m> 


The equations to be solved by Newton-Raphson iteration for the n;. A, , 
and the numbers of moles, are listed in the report by Gordon and McBride. 
Details on initial estimates for the chemical composition, criteria for the 
covergence of the numerical solutions, tests for the possibility of con- 
densed phases, and phase transitions, and mathematical singularities are 
also to be found in the report. 


The thermodynamic data provided in the report with the program listing 
are mostly from the JANAF thermochemical property tables. For each, chemical 
species, the specific heat for constant pressures, the enthalpy, and the 
entropy were computed from the following equations. 
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(35) 


R 


a | + a 2 T + a-^T 2 + a^T- 5 + a^T 


(36) 


£ 

RT 


2 a 3 9 a 4 3 a 5 4 a 6 

aj + — T + — T 2 + — T 3 + — ' T + 

2 3 4 5 T 


(37) 


5 

R 


„ a 3 9 a 4 -z a 5 4 

a .An I + a 2 T + — T t — T 3 + —t + a. 


where the a^ are ’least-squares' coefficients, which are provided in the 
program data deck. Two sets of coefficients are listed, for the tempera- 
ture ranges 300 °K to 1000 °K, and from 1000 °K to 5000 °K. The use of 
either set of coefficients will give identical results at a temperature 
of 1000 °K, 
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DISCUSSION OF RESULTS ed sm 

• ;-!> f\ v ( rp o - j no taeM •, c * A to B'f 
s 4 on 1 wc ! I ct srH' r; l 



The results of the present set of computations for. #(©, decomposition 
of hydrazine are presented in both graphical and tabular form. The tem- 
perature range, from 300 °K to 1 000° K,~i3' divided i nto i ncrements'of "25 ~ 0 Tf. 
The pressure range, from 1.0 atm to 610. atm, is divided Info i rjcfemepts ,of 
10 atm. The tabulated results are given in the appendix'. The graphical f 
presentation consistS-Of three sets of figures, each set consisting. ..of. 
three figures, one, for each of. the react ion product . species,.. The first 
set of figures are plots of mole-f ract io,n (of H 2 ;, or NH^) vs_; temper- j 
ature with pressure as the parameter. The second set of j f igures cons i sts s 
of graphs of mole-fraction vs pressure, with temperature j as {the parameter. 
The fast set of figures consists; of plots of temperature vs pressure with 1 , 
mole-fractions as the parameter.’ ! i i 

For low temperatures and pressures, the decomposition products are 
approximately 2 % .20$ IT?, .and 78$. .Nlfj, _ For pressures greater 'than .5.0 

atm, no H 2 is produced, and the product mixture is 20 % and 80$ HH 3 , . 

For higher temperatures (near 1000 °K) and low pressures, no NH^ is left 
after the decomposition reaction has. finished, and ..the product mixture 
consists of approximately 33$ ^ and 67$, ^ . -For , high temperatures and ....■ 
pressures, the product mixture consists of approximately 59 % H^, 32 % N 2 , 
and 9$ NH 3 . 

Trace amounts of NH and NH 2 may be .found in the gaseous product . . . .. 
mixture for elevated pressures an,d temperatures, but the amounts. of these., 
speci.es are negligible. r '. ; \ 

. ., I t i s cl ear that the ’effect of i ncreas i.ng the pressure, for a given , . 
temperature, is to decrease the amount of F^. Cand subsequently the .amount r 
of N 2 ) in the product mixture, and i ncrease the amount, of NH 3 . .... 


At. low temperatures, the oyer-al l 
represented by 


decomposition 


reaction .may be. 


h « * ^ - y , - < ■ ;r, <• .t , y y * 5 f -yy ' ; 5 y , £ y i : ! ' 3l y , y H Or O'". C ']' v , f, \ A V t c: ■ J '3 L) H • TO 0 f ' f 

At h i gher temperatures ,, the over-a 1 1 reaction may be represented by 

t:nep ; .-iq .y bi-nob y ear brie ;■ .-nutiineqisief ert 'voT tsH; tooqxe 

j- .q&c-.i'; lO'-t'JeQ V • :i i- | t/Tbf cCb.'S bVndad ,yi Ofis <4: i;Q i 7SO j t<?9Vf, i 

oic.-iutfic ■ -. * ■ '<'^2 • : * • ’^2 * no ...HI ,nci tsrutse 


svofl'ei. 

( 2 ) 


The I ow-temperature reaction (equation ( I ), .) is in -complete agreement ; y, 6 
with the results of the i nvestigation made by C. F.> Sayer ( 1 970) . For 
either case, the reaction is exothermic, and proceeds spontaneously. 

Thus , any equilibrium calculation , such as that given here, will show 
no N 2 H 4 concentration for any temperature . 
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A comparison may be made between the results of the present computa- 
tion and the results of A. S. Keston (1967). A summary of the results 
compared is given in the following table. 


Table 5. Comparison of Keston results with present computations. 


P, atm 

T, °K 

Mole Fraction ( 

Keston) 

Mole Fraction (ores 

.ent work) 



h 2 

n 2 

nh 3 

h 2 

n 2 

nh 3 

7.59 

1000 

.55 

(.59) 

.31 

.14 

.665 

.335 

.0012 

32.6 

1054 

.55 

(.56) 

.31 

.14 

.660 

.330 

.0059 

71.0 

1054 

.57 

(.54) 

.32 

,12 
( . 1 5 ) 

.665 

.328 

.01 10 


The numbers in () indicate experimentally-determined values (Keston, 1967) 


A comparison between the sets of results shows good agreement as to 
the amount of nitrogen produced. The discrepancies between the results 
are for the mole-fractions of hydrogen and ammonia. For both the theo- 
retical predictions and experimental findings, Keston' s results i ndi cate 
that fewer moles of hydrogen will be produced, and many more moles of 
ammonia. Qualitatively, both investigations showed that an increase in 
the pressure for a given temperature will result in the production of 
more ammonia. It should be noted, however, that Keston 's results are for 
steady flow of hydrazine across a Shell 405 catalyst bed. The data he 
presents do not appear to have reached zero rate of change with axial posi 
tion in the catalyst bed and thus can be expected to show the lower H 2 
concentrations than would equilibrium calculations; this is because most 
of the H 2 is produced by an endothermic decomposition of NH 3 which pro- 
ceeds at a finite rate, evidently, slower than the transit time of the gas 
through the catalyst bed. 

Two factors, in addition to the basic thermodynamic data used in the 
computation, certainly affect the accuracy of these results: The first is 

the condensation of one or more of the species at certain conditions and 
the second is the real gas behavior of the various species. One would 
expect that for the temperatures and pressures considered in the present 
investigation and N 2 behave substantially as perfect gases, far removed 
from saturation. NH^, on the other hand, may easily saturate (see Table 6). 
Thus one should also consider the non-ideal nature of NH^ for pressures 
and temperatures immediately above saturation in estimating the mole 
fraction of the various species present. 
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Table 6 . Saturation curve of ammonia. 


T 'v °K 

P 'v atm 

195.42 (T.P.) 

0.05997 

200 

0.08553 

210 

0.1756 

220 

0.3347 

230 

0.5982 

240 

1.0124 

250 

1 .6320 

260 

2.526 

270 

3.769 

280 

5.446 

290 

7.652 

300 

10.485 

310 

14.063 

320 

18.49 

330 

23.90 

340 ; 

30.42 

350 

38.19 

360 

47.40 

370 

58.14 

380 

70.60 

390 

84.94 

400 

101.5 

405.6 (C.P.) 

11 1.5 


Figure 10 compares the ideal NH 3 gas pressure for given density and 
temperature with the real gas pressure (Din, 1956) at the same conditions. 
The trends shown in the figure can be explained as follows: As the pres- 

sure increases at constant temperature the energy in the NH 3 molecule 
first tends to migrate from the translational degrees of freedom to the 
rotational or vibrational degree of freedom. Consequently the molecule 
transfers less momentum (exerts a smaller pressure) as the result of elas- 
tic collisions with other molecules than would be expected from a know- 
ledge of its energy content. As one attempts to crowd more and more 
molecules of the same energy content into a given volume, collisions 
finally become so numerous that all the atoms in the molecule become 
involved more or less simultaneously, giving one the appearance of reduc- 
ing the activity in the vibrational degree of freedom. This process is 
most apparent where the total energy content is small, e.g. at low 
temperatures, just above condensation. In the worst case the real gas 
pressure is about a factor of three lower than the ideal gas pressure 
for the same specific volume and temperature. 
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The results cannot be applied exactly to a mixture. The total pres- 
sure gives an indication of the crowding of the molecules but and H2 

molecules may induce a somewhat d liferent response from NH^ than NH^ mol- 
ecules on NHj. General ly, -one would expect that If the NH-* in a mixture 
would suddenly go from a perfect gas to a real gas, the total pressure 
would drop somewhat (in proportion to the change in NH, pressure for the 
prevailing temperature and specific volume) and that the volume fractions 
of N2 and H2 would thereby i ncrease somewhat. However, the total pres- 
sure is then somewhat lower and the NH-, mole-fraction would then tend to 
increase. The net result is probably little different from that given by 
the perfect gas computations. 

Otcd-, i > vC’T 

In comparing equilibrium computations with results for catalytic 
beds it is instructive to '-recall that an adiabatic catalyst can promote 
the redistribution of the energy brought to a reaction by the constitu- 
ents but it does not add or withdraw energy (unless heated or cooled). 

In the case of hydrazine the initial reaction ' 


(3) 3N 2 H 4T *• 4NH 3 + N 2 \ 

is exothermic. If the energy released in this process is all used to 
effect the endothermic decomposition of NH-z to and H 2 then when the 
temperatures and pressures of at I the reaction products equals that of 
the incoming I^H^, the mole fractions of the product species will be as 
given by the equilibrium calculation. Viewed in -this light the equili- 
brium calculations provide an upper 4 bound for the effectiveness of a 
catalytic reactor. ' tv- 


In the real world, of course, things are not so simple. The reactor 
may have to be cooled, reduci ng the amount of heat available to decompose 
the NH^. The reactions may not be complete when the constituents leave the 
reactive area and enter a region of lower pressure and lower temperature." 

In this region the reactions may freeze leaving the mole-fractions what 
they were leaving the reactive area or they may go to completion at the 
lower temperatures 1 and pressures. In the latter ease the NHj concent ra- - 
tions- will be greater and the H2 concentration lower 4 than if the reactions 
froze. Thus to get large quantities of H 2 it is necessary to keep the 
reactants hot long enough for the NHj to decompose completely and then 
cool the mixture rapidly. Another factor to be considered in designing 
a read-tor ¥dr e 'H 2 production from N2H4 is the data of Sayer which indicated 
that a' rhodium catalyst is more effective at promoting direct H 2 formation 
than the iridium (Shell 405) catalyst which seems to favor ammonia pro- 
duction as an intermediate step. 5 5 
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Figure 2. Equilibrium concentraf i on of N~ as a function of temperature at various pressures. 






:entra+ ion 
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Figure 5. Equilibrium concentration of N 2 as a function of pressure at various temp 



pressure at various tei 





MOLE FRACTION 



MOLE FRACTION 



Figure 8. Variation of N 2 pressure with temperature for constant mole fraction 



MOLE FRACTION - NH 



Figure 9. Variation of NH^ pressure with temperature for constant mole fractions 
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perfect gas. 



Comparison of the behavior of NH3 with a perfect gas. 
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APPENDIX A - Nomenclature 


A Helmholtz free energy 

a j j stoichiometric coefficient 

b. ^ a. .n. 

1 ' J J 

b° specified number of kg-atoms of element i per kg of total reactants 
c y constant volume specific heat 
c constant pressure specific heat 

G^ Gibbs free energy function 

g Gibbs free energy per kg of gas mixture 
H enthalpy 

K equilibrium constant, general reaction equation 
k reaction rate constant 

k’ Boltzmann constant 

kj reaction rate constant for first order reaction 

k 2 reaction rate constant for second order reaction 

m number of species in the gas mixture 

N total number of moles 

n number of species in the reacting system 

n^ number of moles of specie X 

ny number of moles of specie Y 

n^y number of moles of reaction product XY 

P pressure 

R gas constant 

S entropy 
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T temperature 

t time 

U internal energy 

V spec i f i c vo 1 ume 

Z partition function 

partition function for specie X 
Zy partition function for specie Y 

Z^y partition function for the reaction product XY 
a thermal expansion coefficient 

3 compressibility 

A Lagrangian multiplier 

y chemical potent i a I 

p density 

Superscripts 

0 quantity at equilibrium 
Subscr i pts 

c value of quantity at the critical point 

1 chemical specie indexing number 

j chemical specie indexing number 

r 'reservoir’ condition 

T value for specified temperature 

CO 3 brackets containing chemical symbols refer to concentrations of those 
■ chemical species. 

Note: Duplicate symbols and/or other symbols appearing in the report are 

defined in the section where used. 


54 



APPENDIX B - Tabulation of Numerical Results 
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CO 

?h 

|h M* 

cn 

00 

CM 

M3 

O' 

© 

O' 

pH 

in © 

© O' 

4 

© 

© CO 

CM 

© 

h* © CO 


O 

o 

h co 

© 

© 

c- 

tn 

in 

fn 

O' 

CM 

CO CO 

CM © 

>4* 

© 

4 

© 

c** 

00 O' © 

CM 

o 

© 

o o 

© 

fH 

HI 

CM 

CO 

* 

in 

?h 

ao O' 

© © 


H 

CM CM 

CM 

CM 

<N nj co 

x 

CM 

CM 

M (NJ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM CM 

CO CO 

CO 

CO 

CO CO 

CO 

CO 

co co tn 


• 

0 

t 0 

0 

• 

• 

t 

• 

0 

« 

t 

* 0 

• 0 

• 

* 

• * 

t 

• 

♦ * 0 


O 

O 

0 

0 

O 

© 

© 

O 

Q 

© 

© 

a 

© © 

© © 

© 

© 

© © 

© 

© 

© © © 


co 

© 

O © 

HI 

CM 

© 

CM 

tn 

© 

in 

CM 

O' so 

?H 4 

© 

4 

© CM 

O' 

CM 

m -fr co 

co 

00 

tn 

O O' 

CM 

in 

© 

M3 

4 

4 

© 

fn 

© © 

f*» -a 

CO 

© 

O' 4 

CM 

4 

N <V O0 

X 

& 

O' 

O' C- 

MO 

co 

O' 

4 

© 

tH 

4" 

© 

O' CO 

© 4 

pH 

© 

a in 

•4* 

CO 

CM CM H 

X 

ih 

fn 

h* fn 

P- 

e- 

M3 

© 

tn 

tn 

4 

CO 

CM CM 

«H HI 

«HI 

© 

© © 

© 

© 

© © © 


• 

# 

* « 

« 

• 

0 

« 

• 

6 

• 

0 

0 • 

• * 

f 

• 

9 * 

* 


0 0 0 


o 

O 

© a 

© 

a 

© 

© 

© 

© 

© 

© 

© © 

0 

0 

© 

© 

© © 

© 

© 

o © © 


4 

C* 

co o 

in 

a 

© 

*H 

tn 

O' 

tn 

© 

© O' 

rsi r*. 

O' 

© 

*h tn 

O' 

HI 

0* O 


H 

pci 

oo h* 

rH 

4 

M3 

© 

O' 

c*» 

O' 

© 

fn fn 

o 4 

pH 

CO 

O' »H 

© 

© 

co » vt 

CM 

O 

O 

© H 

CO 

tn 

© 

CM 

fH 

co 

O' 

© 

pH © 


fn 

O' 

© CM 

CO 

<o 

<h m 

X 

o 

O 

© © 

© 

© 

© 

H 

*H 

CM 

CM 

CO 

4 MC 

in m 

tn 

tn 

© © 

© 

© 

■O 'O «o 


• 

• 

• 1 

* 

• 

• 

• 


« 

• 

• 

* 0 

• « 

• 

• 

• * 

• 

• 

• • • 


o 

o 

© © 

© 

© 

© 

© 

© 

© 

© 

© 

© © 

© © 

© 

a 

© © 

© 

a 

o o o 

a 

o 

in 

© m 

© 

m 

O 

in 

2 

in 

o 

m 

© m 

a m 

© 

m 

© tn 

o 

tn 

o in o 

X 

© 

CM 

ui r- 

© 

cm 

m 

c* 

© 

CM 

tn 

c- 

© CM 

in 

© 

CM 

tn c- 

© 

CM 

iflhO 

UJ 

h 

co 

P0 

co co 

4 

4 

4 

4 

in 

tn 

in 

m 

© © 

© © 

c* 

h* 

P** C"* 

© 

© 

CO CO O' 



PO 

© 

O' O' 

CO 

4 

© 

o 

CM 

© 

tn 

© 

© © 

in r» 

tn 

H 

O' O' 

tn 

c* 

'O co O' 


© 

a 

H CO 

c* 

CM 

O' 

O' 

© 

CM 

m 


O' © 

© CM 

e* 

pH 

co in 

n 

© 

O' o © 

CM 

© 

© 

© © 

© 

H 

fH 

CM 

4 

in 

© 

f*» 

00 O' 

© H 

pH 

CM 

CM CM 

CM 

CM 

CM CO PO 


CM 

CM 

CM CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

nj (Nl 

co ro 

CO 

CO 

CO CO 

co 

PO 

CO PO PO 


0 

0 

0 0 

0 

« 

• 

0 

t 

0 

0 

0 

0 0 

9 0 

0 

9 

# 0 

0 

0 

0 0 0 


© 

© 

© © 

© 

© 

O 

© 

© 

© 

© 

© 

o © 

o o 

© 

© 

0 

0 

© 

© 

O o o 


O 

O' 

’4* *4* 

4 

M3 

P0 

H 

© 

HI 

pH 

O' 

© CO 

fsl O' 

CM 

CO 

eo <o 

O' 

© 

fOHCO 

CO 

© 

4 

© M3 

M3 

in 

CM 

© 

© 

4“ 

?h 

CM 

m © 

rH CO 

in 

CO 

© 4* 

4 


CM 00 

X 

O' 

O' 

00 1- 

tn 

CM 

© 

CM 

m 

© 

o 

CO 

© © 

© CM 

O' 

pH 

m 4 

© 

CM 

CM H H 

X 

c- 

fn 


fn 

c* 

© 

© 

in 

4* 

4* 

CO 

CM CM 

H H 

a 

© 

© o 

© 

© 

© © © 


0 

0 

• § 

0 

0 

0 

0 

# 

0 

0 

0 

0 0 

9 0 

0 

0 

0 0 

0 

0 

0 0 0 


© 

© 

0 

0 

© 

© 

O 

© 

© 

© 

© 

© 

© © 

0 

0 

© 

© 

© o 

© 

a 

o © o 


fH 

CO 

fn M3 

4 

a 

H 

O' 

© 

CM 

4* 

CO 

CM H 

CO 4* 

rO 

© 

CO f- 

© 

tn 

•H SO (O 


pm| 


O' O' 

© 

CM 

© 


?H 

ro 

IH 

O' 

in co 

cm co 

c* 

tn 

O' O' 

r- 

PO 

co^'t 

CM 

a 

© 

O -< 

CO 

M3 

O' 

4* 

Q 

© 

CM 

© 

4* O' 

co © 

© 

© 

H CM 

co 

4 

m- m tn 

X 

© 

© 

O O 

© 

© 

© 

H 

CM 

CM 

CO 

© 

4 4* 

in in 

m 

© 

© © 

© 

© 

>0 lO vO 


0 

0 

t 0 

* 

« 

0 

« 

0 

0 

0 

0 

0 0 

• 0 

0 

0 

0 0 

0 

0 

• * * 


© 

© 

o o 

o 

O 

© 

© 

© 

o 

© 

© 

O Q 

0 

0 

© 

© 

© © 

© 

© 

0 

0 

0 

CL 

© 

tn 

© in 

© 

tn 

© 

tn 

© 

in 

o 

tn 

© tn 

o tn 

© 

tn 

© tn 

© 

tn 

o tn © 

X 

© 

CM 

in h- 

© 

CM 

tn 

t* 

© 

CM 

in 

fn 

© CM 

in 

© 

CM 

in h» 

© 

CM 

me**© 

UJ 

H 

CO 

CO 

CO CO 

-4* 

4 

4- 

4* 

in 

tn 

tn 

in 

© © 

© © 


e* 

fn 

© 

© 

oo ao o> 



4 

© 

fO 1*- 

fH 

© 

HI 

© 

csl 

m 

© 

HI 

© © 

tn 

O' 

© 

© 

©fH 

O' 

© 

mom 


© 

H 

CM 4 

© 

4 

CO 

co 

© 

O' 

CM 

tn 

in 

pH 

© 

© 

ro 

© fH 

© 

O' 

q h »H 

CM 

o 

© 

© © 

© 

tH 

CM 

PO 

4 

tn 

fH 

© 

O' © 

pH 

pH 

CM 

CM 

CM CM 

CM 

CM 


X 

CM 

CM 

CM CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM CO 

CO 

m 

CO 

CO 

CO CO 

CO 

PO 

m m m 


0 

0 

0 0 

0 

0 

0 

0 

0 

0 

0 

t 

0 0 

0 

0 

0 

9 

9 « 

« 

0 

* 0 9 


© 

© 

© Q 

o 

O 

© 

© 

Q 

© 

© 

© 

o o 

© 

o 

© 

© 

© © 

© 

© 

0 

0 

0 


© 

© 

O © 

fH 

't 

*H 

CO 

© 

fH 

© 

CM 

pH fO 

fH 

© 

O' 

pH 

© © 

fH 

fH 

O' r- eg 

CO 

e- 

CO 

© H 

fH 


pH 

fH 

CM 

CM 

CO 

O' 

tn cm 

© 

© 

4 

fH 

fO 

© 

«H 

© CO pH 

X 

O' 

O' 

© fH 

4 

iH 

© 

O' 

CM 

4 

© 

© 

cm r- 

CO 

O' 

fH 

tn 

4 c0 

CM 

CM 

H HH 

X 

c* 


C* 


fH 

© 

m 

tn 

4 

CO 

CM 

CM pH 

H 

© 

© 

© 

© © 

O 

© 

© © © 


0 

0 

• t 

0 

0 

0 

0 

0 

0 

0 

t 

0 0 

0 

0 

9 

t 

0 f 

0 

0 

0 0 0 


© 

© 

o o 

© 

© 

o 

a 

© 

© 

© 

© 

o © 

© 

© 

© 

© 

© © 

© 

© 

© O © 


o 

CM 

e* c* 

© 

© 

fH 

O' 

© 

fH 

CM 

PH 

HI © 

fH 

CO 

CO 

pH 

pH CO 

4 

H 

>o m m 


CM 

m 

pH PO 

CO 

CO 

in 

© 

pH 

fH 

4 

m 

O' co 

fH 

4 

4 

O' 

© © 

4 

O' 

cm in Ph 

CM 

© 

© 

fH CM 

4 

fH 

tH 

© 

CO 

O' 

© 

CM 

fH Csl 

m 

© 

O 

pH 

CO CO 

4 

4 

m m m 

X 

© 

© 

© © 

© 

o 

H 

pH 

CM 

<M 

P0 

4 

4 m 

m 

tn 

© 

© 

© © 

© 

© 

© © © 


0 

# 

0 0 

■ 0 

0 

0 

0 

0 

0 

0 

0 

0 0 

0 

0 

0 

0 

• 0 

0 

0 

0 0 0 


© 

© 

o © 

© 

© 

a 

© 

© 

o 

© 

© 

© © 

Q 

© 

© 

o 

© © 

© 

© 

0 

0 

0 

CL 

© 

m 

© tn 

© 

tn 

© 

tn 

© 

in 

© 

in 

a m 

O 

m 

Q 

m 

© m 

O 

in 

o m o 

X 

o 

CM 

in C 1 * 

© 

CM 

m 

fH 

© 

CM 

m 

fH 

© CM 

in 

fH 

© 

CM 

in fH 

© 

CM 

in ih © 

UJ 

H 

CO 

CO 

co co 

4 

4 

4 

4 

in 

in 

m 

m 

© © 

© 

© 


fH 

fH fH 

© 

© 

ao oo O' 


56 


925 0.6590 0.0092 0.3318 925 0.6565 0.0122 0.3313 925 0.6540 0.0152 0.3308 

950 0.6603 0.0077 0.3320 950 0.6582 0.0102 0.3316 950 0.6561 0.0126 0.3312 

975 0.6613 0.0064 0.3323 975 0.6595 0.0085 0.3319 975 0.6578 0.0106 0.3316 

1000 0.6621 0.0054 0.3324 1000 0.6606 0.0072 0.3321 1000 0.6592 0.0090 0.3318 



OO^M>tCO^O'®N«>OHMHO'^H^^N't«ONQO^O 

fM©©OQO©pHrHCMro4*©r'*©0'0'©rHrH.HfM<MCMiMCM<Mco 

arcMcMCMcMCMCMrMCMCMCMeMCMCMCMcMCMcofOfOfocofOfOcnfOfOfn 

*••**•»*•*#*•*«***•*•#•••*» 

ooooooooooooooooooooooooooo 

* 

ac 

b* r**©©©©f'*<M4 , «*H0',H©O'4*r'*0'fO4 p "H4'CM5n©©h.©0' 

< rO©©rsiinnj*HfM<M<Mfva^©or*QpHrH©cM,HinfMCM4-©foa' 

i^o>^coh^NflomKornNOtf\o ©cm o<0'0^^^n(mh 

O ZhN^NMsh>0^lAinsf^«1NNHP»HOOOOOOOO 

• :#•*•* t g t • ••««*»•**•*##••# # • 

O OOOOOOOOOOOOOOOOOOOOOOOOOOO 

CO 

ONMiA<nmcooa'moN^inh' 4 , Nh^<cmo'H'«hodH 

t\ fM©aOpHCM4‘©0'fo©4’©infHtn0'coin©a'pH<McofO4*4’tn 

tOOOOOOOOHHN(C^^^^intfM^IA^^OO>0'0^ 

a: ooooooooooooooooooooooooooo 

o 

CO 

t/> aO^OWOlAOWO^OinOwOiflOiAOtAOtnOiAO^O 

uj xocMin^ocMinhpocMtnr-^cMtnr-ocMinr^ocMinr^ocMtn 

a: U4mffim<O't«t^'tJn5A^iA4^'O4NNhN»fl0«<O0Na'0' 

a H 


f 


X 

K 

< 

o 

• 

o 

f- 


w 

UJ 


CO 

co 

UJ 

oc 

a 


(M©cnrp#HcotHtMH4tn©©o*4-r , "tn0'©pHCMt i '*©©cM«H©4“ 

ooHfM^eo^^ooHm^ift’tNflo^ooH^^h.Qoo'Cho 

<M©©©©©©pHCM<o4 ,, in©f'-©a'©©pHpH(MCMCM<MCMCMiM<n 

3£rM<MCM<MiM<MCM<NCM<McM<MCMCM<Mio<o<o<nfOfOfomcofoco<o 

i f i « i t i i « • t « i « I i I I f « # « « • I t • 

OOOOOOOOOOOOOOOOOOOOOOOOOOO 


co co >0 cm co O' i** m cm O' n. co co o r* ph in © m *h cm © © © m f« r* 

XC'^0'a5'0 4'HNHlAOO^incom<»'#HC'h*iA^^mNNfH 
X ^ t* fs> r* f* o ^ tn * *0 CM CM r« *4 pH o o o o o o © o o 

*«»*•# *«•*««««« I t « I * * «•#••• 

ooooooooooooooooooooooooooo 


•HO't'*©*©coo'4*4’co©r*cM4*4'inpHfOiPMco©<M©©cMpH 
*H<M©fOinro©to©cMO'^4*? s -cocM4 , ©©©©f , '*tnrMinoviM 
<M©©©*HCM4 , ^©tn©inpHr*<M^rH4*r’*a'©pHlM<0 4’4*4 i tn 

ooooooooooooooooooooooooooo 


aoiAo^OiAOtnoifiOino^OtfiOiAO^oiAOtnoifto 

aEOfNiin^orMUM^o<MuM^OfMtnh-OfMinh»OfNHrvf^oo4in 

H 


X 

H* 

< 

o 

o 

*0 


H 

UJ 

a: 

o 

CO 

CO 

UJ 

oc 

Ou 


CM0inO004 , HN»OininO'NH0»OCy|O'h»O s 0^^N<OfiC 

ooHmincM^mcgtoinhco^ooinH^^N^oNoo^oo 

rvJOOOOOO»<-»csjrO^lft^0f^COOO*^^MfM<NITMM<SI<N«fnm 

xcMCMCMCMcMojcMCMCMCMCMCMCMninifococOforofocococococofO 

• # • * • * i t i « « • « « i • i i * i i i § • i • i 

ooooooooooooooooooooooooooo 


tft*H'H0'fOpH4’CMpHCMpHCM©©fHr*©©©©©inr'''»©©pHfH 

2:hh.h.hNhf%^^iA4 , W^NN^HHO0OOOOOOO 


o 

o 

O 

o 

o 

o 

© 

© 

© © 

© 

© 

© 

©o 

o 

O 

© 

© 

O 

© 

© 

o © 

©© 

© 

co 

CO 

4 * 

«H 

ph 

CM 

CM 


CO CM 

4 * 

co 

m 

CM CO 

CO 

*0 

CM 

in 

© 

tn 

O' 

£ 

in © 

pH 

pH 

CO 

f- 

in 

CO 

co 

h* 

in 

CO © 


r** 

4 - 

© © 

f** 

r- 

H 

O' 

CO 

4 - 

CM 

O' 4 - 

© pH 

mt 

CM o 

o 

O 

pH 

CM 

4 * 


HI 

>0 pH 


ro 

O' 

4 * O' 

CM 

in 

© 

O' 

pH 

CM 

CO 

CO 4 * 

4 * in 

m 

X © 

o 

o 

o 

a 

o 

© 

pH 

pH CM 

CM 

CO 

co 

4 - 4 * 

m 

in 

in 

in 

*0 

*0 

*0 

© © 

© © 

© 

t 

* 

♦ 

• 

• 

• 

• 

• 

« * 

* 

* 

• 

• e 

# 

# 

• 

• 

• 

• 

t 

♦ t 

• • 

• 

© 

o 

o 

o 

o 

© 

© 

© 

o o 

© 

© 

© 

O © 

© 

© 

© 

© 

© 

© 

© 

©Q 

©© 

© 

a, o 

in 

o 

m 

o 

in 

O 

m 

© in 

© 

in 

o 

m © 

in 

Q 

tn 

© 

in 

© 

in 

© in 

o tn 

© 

x o 

<M 

in 

h* 

o 

CM 

in 

T* 

©CM 

m 


© 

cm in 


© 

CM 

tn 

r* 

© 

CM 

in 

© CM 

m 

uj ro 

ro 

CO 

ro 

4 - 

4 - 

4 - 

4 > 

m in 

in 

m 

© 

*OM> 

*D 

r» 

f- 

r- 

r- 

CO 

© 

© © 

O' O' 

O' 


57 


975 0.6561 0.0127 0.3312 975 0.6544 0.0147 0.3309 975 0.6527 0.0168 0.3305 

1000 0.6577 0.0108 0.3315 1000 0.6562 0.0125 0.3312 1000 0.6548 0.0142 0.3310 



PRESSURE * 90.0 ATM. PRESSURE * 100.0 ATM. PRESSURE = 110.0 ATM 



CM 

4 

o 

o 

00 MO 

h- 

CM 

4 

HI 

o 

MO 

M t 

o c* 

in 

O 

CO 

m 

r** 


c- 

O' 

in 

O' 


O 

o 

H 

CM 

co mo 

O 

MO 

co 

CM 

CM 

CM 

m 

<n 

CM 

o 

MO 


in 

O* 

H 

CO 

in 

MO 

CM 

O 

o 

o 

O 

o © 

H 

H 

CM 

CO 

4 

in 

MO 

c- 00 

O' 

o 

O 

HI 

HI 

H 

CM 

CM 

CM 

CM 

X 

CM 

CM 

CM 

CM 

CM CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM CM 

CM 

CO 

CO 

CO 

CO 

CO 

co 

CO 

CO 

CO 


9 

* 

* 

• 

t * 

t 

t 

• 

t 

• 

« 

* 

• « 

9 

9 

9 

9 

« 

9 

9 

9 

• 

9 


O 

o 

O 

o 

o o 

O 

O 

o 

O 

O 

O 

O 

o o 

O 

o 

o 

o 

© 

© 

© 

© 

© 

© 


o 

* 4 * 

o 

00 

cm in 

o 


4 

4 

CM 

m 

4 

04 ITS 

H 

O' 

o 

OO 


m 

in 

© 

|S* 

MO 

m 

O' 


4 * 


r- o 

MO 

CM 

Ch 

c- 

00 

4 

O' 

© C- 

m 

O' 

CM 

o 

in 

n 

O' 

MO 

MO 

© 

X 

O' 

O' 

O' 

00 

f*» sQ 

CO 

O 

in 

o 

4 

00 

H 

m O' 

4 " 

O' 

MO 

fO 

o 

00 

MO 

in 

4 

CO 

X 

c- 

r- 

f- 

c- 

C- 0* 

N 

0* 

© 

MO 

in 

4 

4 * 

CO CM 

CM 

H 

H 

HI 

H 

© 

© 

o 

O 

© 


• 

• 

f 

# 

• 1 

• 

• 

# 

• 

* 

• 

• 

* 9 

« 

• 

* 

« 

9 

9 

9 

9 

9 

e 


o 

o 

O 

© 

© O 

o 

o 

o 

o 

o 

o 

O 

© o 

O 

o 

o 

O 

o 

© 

o 

© 

© 

© 


<0 

CM 

o 

H 

O O' 

ro 

CM 

H 

in 

00 

O' 

H 

00 

4 * 

H 


>0 

M? 

4 


CO 

© 

m 


o 

CM 

in 

o 

O' CM 

CO 

H 

0* 

o 

© 

CM 


O' 00 

CM 

o 

H 

c* 

00 

m 

© 

O' 

IS. 

4 

CM 

o 

O 

o 

HI 

hi ro 

m 

00 

H 

MO 

o 

MO 

*H 

MO hi 

MO 

o 

CO 

m 

C- 

O' 

© 

H I 

CM 

CO 

X 

o 

o 

o 

o 

o o 

Q 

o 

H 

H 

CM 

CM 

CO 


4 

in 

in 

in 

m 

m 

MO 

© 

4) 

MO 


• 

9 

• 

• 

• 1 

9 

* 

9 

* 

• 

• 

♦ 

# • 

9 

• 

« 

» 

9 

9 

9 

« 

9 

9 


o 

o 

o 

o 

oo 

o 

o 

o 

o 

© 

O 

o 

o o 

o 

© 

o 

o 

© 

© 

© 

© 

© 

© 

Qu 

o 

m 

o 

in 

O m 

o 

in 

o 

in 

© 

m 

o 

ir\ o 

m 

o 

in 

o 

n 

o 

in 

© 

in 

o 

X 

o 

CM 

m 


O CM 

in 

c- 

o 

CM 

in 

c- 

o 

n tn 

c* 

Q 

CM 

m 

c- 

© 

CM 

m 

c* 

© 

UJ 

H 

CO 

co 

ro 

ro 

SC 4 

*4* 

4 

in 

in 

m 

in 

4) 

v0 *0 

MO 


h* 

h» 


© 

© 

© 

© 

O' 



CM 

m 

H 

CM 

a © 

CO 

CM 

C* 

r- 

O' 

© 

© 

M* hi 

c- 

© 

H 

© 

o 

CM 

fs. 

© 

CM 

4 


© 

© 

HI 

CM 

4 f* 

H 

r- 

4 

CO 

co 

4 

m 

40 40 

4 

CM 

© 

CO 

c- 

© 

CM 

4 

© 

h- 

CM 

© 

© 

© 

© 

©o 

H 


CM 

CO 

4 

in 

© 

P-* © 

O' 

a 

© 

iH 

H 

CM 

CM 

CM 

CM 

?M 

X 

CM 

<M 

CM 

CM 

CM CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM CM 

CM 

co 

CO 

CO 

CO 

CO 

CO 

CO 

co 

CO 


• 

t 

9 

• 

9 9 

9 

9 

9 

• 

9 

9 

t 

9 9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 


© 

© 

© 

© 

© © 

© 

a 

© 

© 

© 

© 

© 

0 

0 

© 

Q 

© 

a 

© 

© 

© 

© 

© 

© 


O' 

CM 

c* 

© 

© H 

CM 

© 


in 

co 

a 

© 

40 CO 

© 

O' 

4 

c* 

O' 

© 

© 

CM 

© 

CO 

CO 

© 

N 

CO 

fs 

m © 

CM 

40 

H 

c* 

40 

fm{ 

in 

H CO 


fS. 

H 

H 

N 

O' 

4 

CM 

CM 

in 

X 

O' 

O' 

O' 

© 

c» m 

CO 

O' 

in 

O' 

co 


© 

4 © 

fO 

© 

in 

CM 

O' 


© 

in 

4 

co 

X 


C*» 

C* 

r» 


C- 

© 

40 

m 

in 

4 

4 

CO CM 

CM 

HI 

H 


© 

© 

© 

© 

a 

© 


9 

• 

• 

9 

9 9 

9 

9 

9 

9 

9 

9 

9 

9 9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 


o 

© 

o 

© 

0 

0 

© 

© 

© 

© 

© 

© 

CO 

0 

0 

© 

© 

© 

© 

© 

o 

© 

© 

© 

© 


O' 

CO 

CO 

© 

<M O' 

m 

© 

4> 

t** 

h* 

H 

HI 

© 40 

in 

H 

in 

CM 

H 

O' 

fO 

CM 

© 

CM 


o 

CM 

in 

© 

© * 

4) 

40 

CO 

© 

O' 

4 

O' 

CM © 

co 

© 

o 

m 

in 

o 

CO 

CO 

HI 

N 

CM 

© 

Q 

a 

H 

m m 

in 

© 

CM 

40 

•H 

f- 

CM 

© CO 


H 

4 

© 

© 

© 

H 

CM 

CO 

fO 

X 

© 

© 

© 

o 

o © 

© 

© 

H 

H 

CM 

CM 

ro 

CO 4 

4 

in 

m 

m 

m 

© 

© 

© 

© 

© 


9 

9 

9 

9 

# • 

9 

9 

9 

9 

9 

9 

9 

9 » 

9 

9 

9 

9 

« 

9 

9 

9 

9 

9 


© 

© 

© 

© 

o© 

o 

© 

© 

© 

© 

© 

© 

0 

0 

© 

© 

© 

Q 

© 

© 

© 

© 

© 

© 

a. 

© 

m 

© 

in 

o in 

© 

in 

© 

in 

o 

in 

a 

m © 

m 

© 

in 

O 

in 

© 

in 

© 

in 

© 

X 

© 

M 

in 

c* 

O 04 

&n 

c- 

© 

CM 

m 

c* 

© 

cm in 

h* 

© 

CM 

in 


© 

CM 

m 

c- 

© 

as 

t- 

m 

co 

CO 

CO 

-t 

4 

4 

in 

m 

in 

in 

4? 

© © 

© 


C* 

fs. 


© 

© 

© 

© 

O' 



CM 

in 

H 

CO 

rn m 

HI 

CO 

CM 

© 

CM 

4 

in 

H fs» 

pm4 

CM 

© 

© 

4 

CO 

© 


O' 

© 


© 

© 

H 

CM 

4 0* 

CM 

© 

© 

in 

© 

t* 

© 

O' © 

h» 

4 

o 

4 

© 

Hi 

fO 

m 

© 

© 

CM 

© 

o 

O 

© 

© © 

rA 

H 

CM 

co 

4 

in 

© 

c- © 

O' 

© 

H 

H 

H 

CM 

CM 

CM 

CM 

CM 

X 

CM 

CM 

CM 

CM 

04 04 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

04 CM 

CM 

CO 

CO 

CO 

CO 

CO 

ro 

CO 

CO 

CO 


9 

9 

9 

t 

i 9 

9 

9 

9 

9 

9 

9 

9 

9 9 

9 

9 

9 

9 

9 

9 

9 

9 

« 

# 


© 

© 

© 

© 

OO 

Q 

© 

© 

o 

© 

© 

© 

0 

0 

© 

o 

o 

© 

O 

© 

o 

a 

© 

© 


© 

© 

CM 

H 

H CM 

© 

CM 

n* 

-4 

© 

© 

© 

't N 

4 

© 

CM 


© 

CM 

CO 

4 

© 

© 

co 

© 

h- 

CO 

© 

4 m 


© 

CM 

© 

CM 

in 

O' 

in 

h* 

m 

© 

CM 

O' 

CM 

© 


© 

CM 

X 

O' 

O' 

O' 

© 

o- in 

CM 

O' 

4 

© 

CM 

in 

© 

CM © 


N 

4 

H 

© 

h» 

in 

4 

CO 

CO 

X 

r* 

C- 



r** c- 

C- 

© 

© 

m 

m 

4 

CO 

m (\j 

CM 


H 

H 

© 

© 

© 

© 

a 

© 


9 

.9 

9 

9 

9 9 

9 

9 

9 

9 

9 

9 

9 

9 9 

9 

9 

9 

9 

9 

9 


9 

9 

9 


© 

© 

o 

© 

a © 

© 

© 

© 

© 

© 

© 

© 

o o 

© 

© 

© 

a 

© 

© 

a 

© 

o 

© 


© 

in 

f- 

© 

© 4 

CO 

in 

HI 

CM 

© 

© 

in 

in © 

in 

© 

© 

CM 

© 

n 

© 

fH 

CO 

o 


i—4 

CM 

n 

H 

H C- 

© 

.Hi 

H 

© 

sM 

© 

CM 

m ?n 


© 

O' 

CO 

i— I 

© 

© 

h* 

4 

© 

CM 

Q 

© 

a 

H 

CM CO 

© 

O' 

CO 

C* 

CO 

© 

4 

O' 4 

© 

CM 

4 

fs. 

O' 

© 

H 

CM 

CO 

4 

X 

O 

© 

o 

o 

© © 

© 

o 

HI 

H 

CM 

CM 

ro 

CO xfr 

4 

m 

tn 

in 

m 

© 

© 

© 

© 

© 


9 

• 

9 

• 

• 9 

9 

9 

9 

9 

9 

9 

9 

9 9 

9 

9 

9 

9 

9 

9 

9 

« 

9 

• 


© 

© 

© 

© 

© © 

© 

© 

o 

a 

© 

© 

o 

© o 

© 

Q 

a 

© 

© 

© 

a 

O 

© 

© 

a 

o 

m 

© 

m 

© in 

© 

in 

© 

in 

o 

in 

o 

in © 

in 

O 

in 

o 

in 

a 

m 

a 

m 

O 

X 

o 

CM 

in 

N» 

© CM 

in 

c* 

© 

CM 

in 

h* 

o 

cm in 

r- 

© 

CM 

in 

h- 

o 

CM 

in 

r- 

© 

LU 

H* 

CO 

ro 

co 

CO 

4 4 

4 

4 

in 

m 

in 

n 

© 

© © 

© 


N 


fs. 

© 

© 

© 

© 

O' 


58 
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975 0.6460 0*0247 0.3292 975 0.6444 0.0267 0.3289 975 0.6428 0.0286 0.3286 

1000 0.6491 0.0211 0.3298 1000 0.6477 0.0227 0.3295 1000 0.6463 0.0244 0.3293 
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PRESSURE * 210*0 ATM. PRESSURE * 220.0 ATM. PRESSURE = 230.0 ATM 
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PRESSURE * 270.0 ATM. PRESSURE * 280.0 ATM. PRESSURE = 290.0 ATM 
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900 0.5957 0.0852 0.3191 900 0.5935 0.0878 0.3187 900 0.5914 0.0904 0.3183 

925 0.6066 0.0721 0.3213 925 0.6047 0.0744 0.3209 925 0.6028 0.0766 0.3206 

950 0.6156 0.0613 0.3231 950 0.6139 0.0633 0.3228 950 0.6123 0.0653 0.3225 

975 0.6230 0.0524 0.3246 975 0.6216 0.0541 0.3243 975 0.6202 0.0558 0.3240 

1000 0.6292 0.0449 0.3258 1000 0.6280 0.0464 0.3256 1000 0.6267 0.0479 0.3253 
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925 0.5955 0.0854 0.3191 925 0.5937 0.0876 0.3187 925 0.5919 0.0897 0.3184 

950 0.6059 0.0729 0.3212 950 0.6043 0.0748 0.3209 950 0.6027 0.0767 0.3205 

975 0.6145 0.0625 0.3229 975 0.6132 0.0642 0.3226 975 0.6118 0.0658 0.3224 

1000 0.6218 0.0538 0.3244 1000 0.6206 0.0553 0.3241 1000 0.6194 0.0567 0.3239 
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975 0.6105 0*0674 0.3221 975 0.6091 0.0691 0.3218 975 0.6078 0.0707 0.3215 

1000 0.6182 0.0582 0.3236 1000 0.6170 0.0596 0.3234 1000 0.6158 0.0610 0.3232 
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PRESSURE * 450.0 ATM. PRESSURE * 460.0 ATM. PRESSURE * 470.0 ATM 
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PRESSURE * 510.0 ATM. PRESSURE * 520.0 ATM. PRESSURE * 530.0 ATM 
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